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27 Abstract 


28 Measurement capabilities of five meteor radars are assessed and compared to determine 29 how well radars 
having different transmitted power and antenna configurations perform in 30 defining mean winds, tidal amplitudes, 
and gravity wave (GW) momentum fluxes. The five 31 radars include two new-generation meteor radars on Tierra 

o o 

del Fuego, Argentina (53.8 S) and on 32 King George Island in the Antarctic (62.1 S) and conventional meteor radars 
at Socorro, New 33 Mexico (34.1 N, 106.9 W), Bear Lake Observatory, Utah (-41.9 N, 1 1 1.4 W), and Yellowknife, 

o o 

34 Canada (62.5 N, 1 14.3 W). Our assessment employs observed meteor distributions for June of 35 2009, 2010, or 
2011 for each radar and a set of seven test motion fields including various 36 superpositions of mean winds, constant 
diurnal tides, constant and variable semidiurnal tides, 37 and superposed GWs having various amplitudes, scales, 
periods, directions of propagation, 38 momentum fluxes, and intermittencies. 39 Radars having higher power and/or 
antenna patterns yielding higher meteor counts at 40 small zenith angles perform well in defining monthly and daily 
mean winds, tidal amplitudes, 41 and GW momentum fluxes, though with expected larger uncertainties in the daily 
estimates. 42 Conventional radars having lower power and a single transmitting antenna are able to describe 43 
monthly mean winds and tidal amplitudes reasonably well, especially at altitudes having the 44 highest meteor 
counts. They also provide qualitative estimates of GW momentum fluxes at the 45 altitudes having the highest 
meteor counts; however, these estimates are subject to uncertainties 46 of -20 to 50% and uncertainties rapidly 
become excessive at higher and lower altitudes. 47 Estimates of all quantities degrade somewhat for more complex 


motion fields. 48 49 50 1 



Introduction 


5 1 Gravity waves (GWs) play significant roles in atmospheric dynamics, chemistry, 52 microphysics, and plasma 
processes extending from Earth’s surface into the thermosphere and 53 ionosphere. At lower altitudes, GWs 
contribute to boundary layer and mesoscale structures, 54 organize convection, induce turbulence and mixing of 
relevance to weather prediction and 55 aircraft at flight altitudes. In the middle stratosphere, GWs induce polar 
stratospheric clouds that 56 contribute to ozone destruction. Throughout the atmosphere, but especially at higher 
altitudes, 57 GWs systematically influence the large-scale circulation and thermal structure via energy and 58 
momentum transport and deposition. These influences are typically slow and systematic at lower 59 altitudes, but 
may be local and strong accompanying severe events. In the mesosphere and lower 60 thermosphere (MLT), GWs 
lead to closure of the mesospheric jets, large departures from 61 radiative equilibrium, strong turbulence and mixing 
of heat, momentum, and constituents, and 62 significant interactions with, and influences on, tides and planetary 
waves. GWs are expected to 63 have similar, but still largely unknown, effects higher in the thermosphere and 
ionosphere. 64 Indeed, the vertical transport of horizontal momentum by GWs, and its deposition accompanying 65 
GW dissipation, account for the major influences of these motions on the mean and large-scale 66 atmospheric 
circulation, structure, and variability. Reviews of these various dynamics 67 emphasizing the higher altitudes are 
provided by Hines [1960], Lighthill [1978], McIntyre 68 [1989], Hocke and Schlegel [1996], Nappo [2002], Fritts 
and Alexander [2003], Kim et al. 69 [2003], Fritts et al. [2006a], and Fritts and Lund [201 1], 70 The various 
responses to momentum transport by GWs noted above make this a key 7 1 quantity in understanding and modeling 
atmospheric structure at all altitudes. Consequently, 72 significant efforts have addressed these dynamics employing 


theoretical, modeling, and 



73 observational methods, typically focusing on the higher altitudes where the effects of these 


74 dynamics are most pronounced [Fritts and Alexander, 2003], Further quantification of these 75 dynamics 
throughout the atmosphere also remain a major research need, given the important of 76 their proper 
parameterization for weather and climate forecasting [Kim et al., 2003]. 77 Prior to new meteor radar capabilities, 
direct measurements of GW momentum fluxes 78 have only been possible with Doppler radars or lidars having 
symmetric, relatively narrow, of£79 zenith coplanar beam pairs [Vincent and Reid, 1983; Reid and Vincent 1987; 
Fritts and Vincent, 80 1987; Fukao et al., 1988; Reid et al., 1988; Fritts and Yuan, 1989; Fritts et al., 1990, 1992, 81 
2006b; Fritts and Janches, 2008; Sato, 1990, 1993, 1994; Tsuda et al., 1990; Wang and Fritts, 82 1990, 1991; 
Hitchman et al., 1992; Nakamura et al., 1 993; Murayama et al., 1994; Murphy and 83 Vincent, 1993, 1 998; Acott et 
al., 2009] or via in situ measurements by aircraft at lower altitudes 84 or chaff measurements in the MLT [e.g., Lilly 
and Kennedy, 1973; Lilly et al., 1982 ', Brown, 85 1983 ', Meyer et al., 1 989; Nastrom and Fritts, 1992; Smith et al., 
2008], Indirect methods have 86 nevertheless contributed to quantification of, or constraints on, momentum fluxes 
employing 87 airglow [e.g., Swenson et al., 1999; Fritts et al., 2002; Espy et al., 2004, 2006], balloon [e.g., 88 
Hertzog and Vial, 2001; Hertzog et al., 2008], and satellite measurements [e.g., Em et al., 2004], 89 These various 
measurements have revealed a variety of responses, including systematic seasonal90 mean momentum fluxes 
accounting for large-scale forcing at lower and higher altitudes, 91 responses to specific sources and source regions, 
episodic and strong forcing accompanying 92 GWs that are transient and localized, but achieve large amplitudes, and 
significant filtering and 93 interactions with larger-scale mean, tidal, and planetary wave (PW) motions. 94 Ideally, 
routine momentum flux measurements would provide continuous sensitivity to 95 the magnitudes that are 


dynamically important and also capture the temporal variations 



96 accompanying strong, but transient, events that may nevertheless contribute significantly to 
97 hourly, daily, monthly, and seasonal means. Improved understanding and parameterization of 98 GW 
influences also impose a need for such measurements at a range of locations spanning 99 representative source 
regions for GWs and the environments in which they propagate. 100 Unfortunately, those systems that do (or 
did) have such measurement capabilities are very 101 limited, their costs are very high (typically ~$1 to 10M, 
with several much more costly), very 102 few measure continuously, and only two (to our knowledge), sodium 

o 

resonance lidars at the 103 Arctic Lidar Observatory for Middle Atmosphere Research (ALOMAR) at 69.3 N in 

o 

Norway 104 and the Andes Lidar Observatory (ALO) at 30 S in Chile, are even at the edges of regions of 105 
identified high seasonal GW activity extending into the stratosphere and MLT. So there are 106 considerable 
motivations for exploring other measurement techniques that may offer the 107 potential to significantly expand 
such measurements with much cheaper systems. Such systems 108 should also have the potential for continuous 
measurements; optical systems, whether active 109 (lidar) or passive (airglow, FPI, etc.) cannot provide full 
diurnal and annual measurements, 110 independent of weather, suggesting that only radar measurements are 
likely to address this need. 1 1 1 The Southern Argentina Agile MEteor Radar (SAAMER) [Fritts et al., 2010a, 
hereafter 112 FlOa ] was specifically designed to add a potential for GW momentum flux measurements to the 
113 conventional meteor radar capabilities for mean, tidal, and PW measurements for the reasons 114 noted 
above. This was accomplished by devising a radar beam pattern that yielded a majority of 1 15 meteors at zenith 
angles <50 and employing high power to achieve as high a meteor rate as 1 16 possible at these angles. 

o o 

SAAMER was installed at Rio Grande on Tierra del Fuego, Argentina 1 17 (53.8 S, 67.7 W) and became 
operational in May 2008. SAAMER momentum flux measurement 118 capabilities were evaluated using 


observed meteor distributions for September 2008, a series of 



119 test motion fields including prescribed mean winds, diurnal and semidiurnal tide winds, and 
120 GWs having various spatial and temporal scales [Fritts et al., 2010b, hereafter FlOb], and the 121 analysis 
procedure suggested by Hocking [2005] employing a “full-field” fit to the various 122 winds, variances, and 
covariances characterizing the superposed mean winds, tides, and GWs. 123 These suggested that for real wind fields 
having similar composition and character to the test 124 fields, SAAMER should be expected to provide very good 
definition of the mean and tidal winds 125 and reasonably good definition of GW momentum fluxes where meteor 
count rates are sufficient. 126 The apparent success of SAAMER in measuring GW momentum fluxes in addition 
to 127 mean, tidal, and PW winds was the motivation for a SAAMER clone, the Drake Antarctic Agile 128 MEteor 

o o 

Radar (DrAAMER) that was installed at the Brazilian Antarctic Comandante Ferraz 129 Base (62. 1 S, 58.7 W) in 
March 2010. Comparisons of the mean and tidal winds and the GW 130 momentum fluxes measured by SAAMER 
and DrAAMER by Fritts et al. [2011, hereafter FI 1] 131 demonstrated close correspondence between mean and 
tidal winds at the two sites, as well as 132 significant similarities between the inferred GW momentum fluxes 
between years and sites from 133 April to June of 2010 and 2011. This agreement suggested significant confidence 
in the GW 134 momentum flux measurement potential for both radars, given the tests performed with SAAMER 135 
data. Nevertheless, we believe it is important to subject DrAAMER to the same tests applied to 136 SAAMER. We 
also believe that applications of these tests to additional meteor radars having 137 different characteristics 
(particularly power and beam configuration) may help us identify what 138 radar capabilities are required to achieve 
GW momentum flux measurements that are sufficiently 139 accurate to be valuable to the community. For example, 
the Hocking [2005] analysis method was 140 tested with meteor radar data from Socorro, NM and Resolute Bay, 
Canada and judged to 141 provide reasonable two-month estimates (though without validation). Other meteor radars 


have 



142 also been employed for GW variance, momentum flux, and tidal modulation studies and suggest 
143 that more radars may provide enhanced measurement capabilities with testing and suitable 144 analysis 
techniques [Kumar et al., 2007; Antonita et ah, 2008; Clemesha et ah, 2009; Mitchell 145 and Beldon, 2009; Placke 
et ah, 201 1]. A related study by Vincent et al. [2010] employed a 146 Monte Carlo analysis of measurement 
capabilities assuming a conventional meteor radar beam 147 pattern and expected radial velocity and angle-of-arrival 
uncertainties and concluded that 148 accurate momentum flux estimates with such systems would likely require long 
averages. 149 Our goals in this paper are to evaluate SAAMER and DrAAMER GW momentum flux 150 
measurement capabilities relative to several other conventional meteor radars in order to 1 5 1 determine 1 ) what radar 
parameters are required to enable such measurements, 2) whether other 152 radars can also provide credible GW 
momentum flux measurements in their current 153 configurations, or 3) whether system upgrades are needed to 
enable these capabilities. The other 154 radars employed for this assessment are at Socorro (SRO), New Mexico 
(34.1 N, 106.9 W), Bear 155 Lake Observatory (BL), Utah (-41.9 N, 1 1 1.4 W), and Yellowknife (YKF), Canada 

o o 

(62.5 N, 156 1 14.3 W) (see Hocking , 2005, and Hocking and Kishore Kumar , 2011, respectively, for more 157 
detailed descriptions of the SRO and YKF radars). To address these goals, our paper is structured 158 as follows. 
Section 2 provides a summary of the characteristics of the radars employed for this 159 assessment and the analysis 
methods used to evaluate the performance of each. Test motion 160 fields, which include mean winds, diurnal and 
semidiurnal tides, and/or GWs having various 161 characteristics are described in Section 3. The performance of the 
five radars for the various tests 162 employing real spatial and temporal sampling for each is described in Section 4. 
Section 5 163 summarizes our findings and discusses measurement accuracies relative to other assessments. 164 Our 


conclusions are presented in Section 6. 



165 2 


Characteristics of Radars Employed for this Study 

166 We employ five different meteor radars for our assessment of GW momentum flux 167 measurement 
capabilities in this study. Two of these were specifically designed intending to 168 provide this capability (i.e., 
SAAMER and DrAAMER); the other three are conventional meteor 169 radars having lower peak power and a 
single transmitting antenna that had not anticipated such 170 measurements (i.e., the SKiYMET meteor radars 
previously installed at Socorro, NM, Bear Lake 171 Observatory, UT, and Yellowknife, Canada; see Table 1). The 
relative locations of the five 172 radars and the meteor distributions for each are shown for one day during the month 
in Figure 1. 173 Meteor distributions in altitude, by day throughout the month considered, and by hour 174 
throughout a composite day for each radar are shown for comparison in Figure 2. The meteor 175 distributions in 
Figure 1 all exhibit nulls in the polar diagrams that arise from excluding meteors 176 at range increments 
corresponding to the pulse repetition frequency (PRF) in order to exclude 177 contamination due to ground clutter. 
SAAMER and DrAAMER use a smaller PRF, and thus 178 have these nulls at larger spacing. The statistics in Figure 

o 

2 include meteors at zenith angles from 179 15 to 60 for all radars. As expected, meteor counts are significantly 
larger for the radars having 180 higher power, and higher at smaller zenith angles for SAAMER and DrAAMER than 
for the 181 radars with single-antenna TX systems. They also vary significantly throughout the month, 182 especially 
at SAAMER and DrAAMER, and exhibit large diurnal variations at all radars. 183 The GW momentum flux 
measurement potential of SAAMER was previously evaluated 184 for a number of test fields by FI Ob. Comparisons 
of SAAMER and DrAAMER measurements 185 for April, May, and June of 20 1 0 and 2011 were presented by FI 1, 
and suggested that 186 DrAAMER momentum flux estimates provided a consistent picture of GW momentum fluxes 
187 accompanying strong GW sources over the Drake Passage, but did not evaluate the DrAAMER 



188 momentum flux measurement capability directly. An initial application of the method employed 
189 here was also used by Hocking [2005] to estimate two-monthly GW momentum fluxes by 190 SkiYMET meteor 
radars at Socorro, NM and Resolute Bay, Canada. These estimates were seen 191 to be roughly consistent with 
expected values, but were not evaluated in detail. However, a 192 systematic evaluation of the relative momentum 
flux measurement potential of these various 193 radars with specific test fields has not been performed. This is the 
goal here, and the 194 characteristics of the five radars are described for comparison in Table 1. We employ meteor 
195 distributions obtained during June 201 1 for KGI and TdF, June 2010 for SRO and BL, and June 196 2009 for 
YKF, for which meteor counts tend to be larger at northern than at southern latitudes. 197 3 Specification of Test 
Motion Fields 198 We showed in FI 0b that the SAAMER beam pattern and meteor counts enable relatively 199 
high-precision measurements of mean winds and tides over fairly short intervals, and that GW 200 momentum fluxes 
can be estimated with reasonable accuracies where meteor counts are high, 201 even when the large- and small-scale 
motion field is variable on multiple time scales. The first 202 measurements over DrAAMER and comparisons with 
those over SAAMER by FI l suggest that 203 the same can likely be said for DrAAMER. Here we repeat the tests 
previously applied to 204 SAAMER for DrAAMER and the three conventional meteor radars at Socorro, MN, Bear 
Lake 205 Observatory, UT, and Yellowknife, Canada in order to evaluate DrAAMER capabilities more 206 
completely and determine whether, and under what conditions, useful GW momentum flux 207 measurements may 
also be possible with the other convenional meteor radars. 208 Following FI 0b, we employ real meteor spatial and 
temporal distributions observed by 209 SAAMER, DrAAMER, and the three other meteor radars during various 
months to evaluate the 210 measurement capabilities of each. Radial velocities at each meteor location and time for 


each 



2 1 1 radar and each month assessed are specified by each of seven test velocity fields. Each test field 


212 includes superposed mean, tidal, and/or GW velocity fields with constant and/or spatially- and 213 
temporally-variable tidal and GW amplitudes that are intended to be representative of the scales, 214 amplitudes, and 
momentum fluxes in the MLT over SAAMER and DrAAMER. Mean, tidal, and 215 GW parameters defining these 
motion fields for each case are listed in Table 2. The test fields 216 are those employed by FI Ob, range from highly 
idealized and stationary to spatially and 217 temporally modulated at large and small scales, and have the following 
components: 218 U(x,y,z,t) = Um + Uo(z,t) sin(2;rt/TD) + UsD(z,t) sin(27tt/TsD) 219 + UGwi(x,y,z,t) sin(kix + liy + 
nnz - 27tt/TGwi) 220 + Uow 2 (x,y,z,t) sin(k 2 x + by + rmz - 2?tt/TGW2) 22 1 + Uow3(x,y,z,t) sin(k3x + rmz) (1) 222 
V(x.y.z,t) = Vm - VD(z,t) cos(27tt/To) - VsD(z,t) cos(27tt/TsD) 223 + VGwi(x,y,z,t) sin(kix + liy + rmz - 2jrt/TGwi) 224 
+ VGW 2 (x,y,z,t) sin(k 2 x + by + rmz - 2jrt/TGW2) 225 + VGW4(x,y,z,t) sin(l4y + rmz) (2) 226 W(x,y,z,t) = 
W(,wi(x,y,z,t) sin(kix + nnz - 2jtt/TGwi) 227 + WGW 2 (x,y,z,t) sin(k 2 x + by + m 2 z - 27tt/TGW2) 228 + WGW3(x,y,z,t) 


sin(k3x + m3z) 229 + Wow4(x,y,z,t) sin(l4y + rmz) (3) 



230 Note that we employ the same test cases defined by FlOb, but these Eqs. (1-3) correct several 


23 1 typos appearing in the equations presented in FI Ob. 232 Each test field includes some or all of the following 
components: 233 a. zonal and meridional mean winds, Um and Vm, 234 b. diurnal and semidiurnal tides having zonal 
and meridional amplitudes of (Ud, Vd) 235 and (Usd, Vsd) assumed to rotate counterclockwise with time (assuming 
a southern 236 hemisphere location), and which may have either constant or varying amplitudes with 237 increasing 
altitude, 238 c. traveling GWs having amplitudes (Ugw, Vgw, Wgw), zonal, meridional, or oblique 239 
propagation, spatial and temporal variability, correlated horizontal and vertical 240 motions, and constant or variable 
momentum fluxes, and 241 d. stationary mountain waves (MWs) having zonal and meridional propagation, only 242 
spatial variability, correlated horizontal and vertical motions, and constant 243 momentum fluxes. 244 In each case, 
the test field amplitudes were chosen to correspond roughly to measured 245 values over SAAMER and DrAAMER 
on TdF and KGI, which include large semidiurnal tides 246 and GW momentum fluxes [FlOa, FI Ob, Fll], Large 
semidiurnal tide amplitudes at these sites 247 are consistent with expectations of the most recent version of the 
Global-Scale Wave Model 248 (GSWM-09) [Zhang et al., 2010a, b, Fll], Large GW momentum fluxes are 
suggested by the 249 major global hotspot of GW activity in the stratosphere and lower mesosphere centered over 
the 250 Southern Andes, Drake Passage, and Antarctic Peninsula [FlOa, and references therein], 25 1 4 Evaluation 


of Radar Measurement Capabilities 



252 We describe here the performance of each of the five radars for the seven test cases 253 described above. Mean 
winds, diurnal and semidiurnal tidal amplitudes, and GWs having 254 prescribed spatial and/or temporal variability 
and momentum fluxes propagating zonally, 255 meridionally, or at other azimuths are defined by Eqs. (1-3) for the 
seven cases listed in Table 2. 256 Specified and recovered (i.e., the velocity fields inferred from our S-transform tidal 
fits and the 257 Hocking statistical analysis) profiles of the mean winds, diurnal and semidiurnal tides, and GWs 258 
momentum fluxes, and their daily and composite-day hourly variability, where appropriate, are 259 discussed 
separately for each case below. 260 b.l. Case 1 261 Specified fields, and those recovered using observed meteor 
distributions in space and 262 time for each radar for June 2009, 2010, or 201 1 as described above, are shown as 
monthly 263 means for Case 1 in Figure 3. In this case, tides have only temporal variations and GWs have 264 only 
horizontal variations, with 50 and 100 km wavelengths in the zonal and meridional motions, 265 respectively (see 
Table 2). 266 As seen in Figure 3, agreement between specified and recovered fields is highly variable 267 among 
the five radars. As might be expected based on meteor counts, mean winds and diurnal 268 and semidiurnal tide 
amplitudes are recovered quite well at KGI, TdF, and BL. Mean winds and 269 diurnal tide amplitudes are recovered 
within a few percent at each of these sites, with slightly 270 greater uncertainties at BL than at the other two. 
Semidiurnal tide amplitudes are systematically 271 less than specified at all three sites, but only by ~2-3% at TdF 
and BL, and by ~5% at KGI. 272 Mean winds and tidal amplitudes exhibit similar tendencies and accuracies at SRO 
and YKF at 273 the central altitudes, but increasing errors at the lowest and highest altitudes where meteor counts 


274 are small. 



275 GW momentum flux estimates at KGI, TdF, and SRO are likewise accurate within — 10%, 276 except at KGI and 


TdF below 79.5 km. In contrast, SRO and YKF exhibit errors of -20% or 277 larger at intermediate altitudes, again 
with larger errors where meteor counts are small. 278 b.2. Case 2 279 Case 2 includes no mean winds or tides and 
only a zonal stationary GW with prescribed 280 horizontal and vertical phase variations with wavelengths of 50 and 

15 km (not physical, but 281 simple, see Table 2). Referring to Figure 4, we see that, as in Case 1, the recovered 

-1 

mean and 282 tidal motions are quite accurate at KGI, TdF, and BL, typically within -2 ms or better, except 283 for 

-1 

the tidal amplitudes at BL above -92 km. Estimates at SRO and YKF are slightly less 284 accurate, within -5 ms , 
but nevertheless quite reasonable, except at the highest and lowest 285 altitudes, especially at YKF. 286 Unlike Case 

1, all the radars except YKF provided accurate estimates of GW momentum 287 fluxes in this case, with the largest 

2 -2 

departures of only -5 m s at the highest and lowest altitudes 288 shown. However, meteor counts at YKF were 
apparently too small to yield adequate statistics, 289 except at the central altitudes having the highest meteor counts, 
from 82.5 to 88.5 km. The 290 significant differences in the accuracies of the GW momentum fluxes at all sites and 
in the mean 291 wind and tidal amplitudes between Cases 1 and 2 suggest a strong influence of large and variable 
292 winds on the quantification of GW momentum fluxes, even for persistent and well-defined GWs. 293 These 
influences will be quantified further below and discussed in some detail in Section 5. 294 b.3. Case 3 295 Case 3 

considers a motion field varying from that of Case 2 only in having two 296 orthogonal GWs, each of which has no 

phase variation in altitude or time and with the GW 297 propagating in the zonal plane having a large vertical 

-1 

velocity of 10 ms (again not physical, but 



298 simple). Results displayed in Figure 5 show that mean motions, tidal amplitudes, and GW 


299 momentum fluxes are again described well with the meteor distributions for KGI, TdF, and BL, 300 and the 
mean motions and GW momentum fluxes are described reasonably at SRO and at YKF, 301 but with increasing 
errors in momentum fluxes at the higher altitudes. However, tidal amplitudes 302 now exhibit larger errors at SRO 
and YKF than seen in Case 2, with these motions reasonably 303 defined at ~85 to 92 km at YKF and in the zonal 
component at SRO, but poorly defined at lower 304 and higher altitudes at both these sites and even at central 
altitudes in the meridional components 305 at SRO. This cannot be a result of insufficient sampling of the mean and 
tidal motions, as they 306 are all zero. Instead, the cause appears to be the larger GW vertical velocities that are not 
307 sufficiently averaged in defining large-scale (horizontal) winds, despite the validity of the 308 momentum flux 
estimates for stationary GWs at these sites (see further discussion in Section 5). 309 b.4. Case 4 310 We now 
consider a more complex superposition of mean, tidal, and GW fields given by 3 1 1 the sum of the mean and tidal 
motions in Case 1 and both stationary and propagating GWs in 312 both the zonal and meridional planes (see Table 
2). This case differs from the previous cases in 313 that both zonal and meridional GW momentum fluxes now have 
contributions from stationary 314 and propagating GWs having different spatial structures (with horizontal 
wavelengths of 30, 40, 315 50, and 100 km). Case 4 also includes both large semidiurnal tide amplitudes and GW 
vertical 316 velocities that were suggested above to contribute to errors in the estimates of mean and tidal 317 

amplitudes and GW momentum fluxes. Depending on the phases of the four GWs, vertical 318 motions may be ~15 

-1 

ms or larger in this case. 319 Results displayed in Figure 6 show that errors in estimates of the mean winds and tidal 
320 amplitudes are comparable (and small) for Cases 1 and 4 at KGI, TdF, and BL. Semidiurnal tide 



321 amplitudes are very slightly underestimated, with the largest errors of~5% at KGI (note the 
322 different scales in the third columns of Figures 3 and 6). As seen in Case 1, mean wind and tidal 323 amplitude 
estimates at SRO and YKF for Case 4 are also reasonable at central altitudes, but 324 exhibit greater errors at lower 
and/or higher altitudes and at central altitudes in the meridional 325 diurnal tide amplitudes. 326 Turning now to 

GW momentum flux estimates for Case 4, we see that these are very 327 good at KGI, TdF, and BL, with 

2 -2 

comparable errors (~a few ms) except at the highest and 328 lowest altitudes where they are somewhat larger. 

2 -2 

Similar estimates are also obtained at SRO and 329 YKF at central altitudes (with errors for Case 4 of ~10 m s or 
less), with much larger errors 330 occurring at lower and higher altitudes. Apparently, mean winds and large tidal 
amplitudes do 331 not prevent accurate estimates of GW momentum fluxes when meteor counts are sufficiently 332 
high and GWs are sufficiently sustained (whether stationary or propagating) to enable sampling 333 of all phases 
throughout the month. Large GW amplitudes (especially vertical velocities), 334 however, do cause errors in 
estimates of means winds and tidal amplitudes where meteor counts 335 are not sufficiently high. 336 b.5. Case 5 
337 Case 5 explores a superposition of larger mean winds, a larger, but constant, diurnal tide, 338 a semidiurnal tide 
that exhibits both a 10-day amplitude modulation and amplitude growth with 339 altitude, and four superposed GWs. 
The latter include 1) stationary GWs having zonal and 340 meridional orientations, constant amplitudes and 
momentum fluxes, and horizontal wavelengths 341 of 30 and 40 km and 2) propagating GWs having zonal and 
meridional orientations, horizontal 342 wavelengths of 50 and 100 km, respectively, and amplitude (and momentum 


flux) modulation 343 by, and anti-correlations with, the semidiurnal tide amplitudes. 



Results for Case 5 are displayed in four formats. Figure 7 shows monthly mean profiles 
345 of mean winds, tidal amplitudes, and GW momentum fluxes as discussed for Cases 1 to 4 above. 346 Figure 8 
displays daily estimates of mean winds and diurnal tide amplitudes throughout the 347 month in the presence of 
variable semidiurnal tide and GW amplitudes. Figure 9 compares daily 348 estimates of semidiurnal tide amplitudes 
with the specified variations. Hourly estimates of GW 349 momentum fluxes throughout the composite day are 
compared with the specified variations in 350 Figure 10. 35 1 Results shown in Figure 7 reveal mean wind estimates 

■4 

(first column) that are very 352 accurate at all radars (within ~1 ms except at the lowest and highest altitudes at SRO 
and YKF). 353 Diurnal tide amplitude estimates (second column) are likewise very good at TdF, nearly as good 354 

-t . -1 

at KGI and BL (within ~1 to 2 ms ), and reasonable at central altitudes at YKF (within ~2 ms ). 355 Somewhat 

-1 

larger errors (~5 ms ) are again seen in the meridional diurnal amplitude at SRO. 356 Mean semidiurnal tide 
amplitude estimates (third column) are seen to be very precise at TdF and 357 BL, to exhibit slight underestimates at 

KGI, and to have comparable errors about the specified 358 mean amplitudes at SRO. Larger errors are seen at YKF 

-1 

(as large as -5-10 ms ), but the 359 amplitude growth with altitude is still captured at the central altitudes. 360 

Daily mean wind and diurnal tide amplitude estimates shown for Case 5 in Figure 8 are 361 most accurate at TdF and 

-1 

BL (with RMS uncertainties of ~1 ms or less) and somewhat less 362 accurate at KGI (with RMS uncertainties of 

-1 

-2 ms ), except at the lowest altitude at KGI and the 363 lowest and highest altitudes at BL where meteor counts are 
small. Daily mean estimates are 364 considerably less accurate at SRO and YKF due to the significantly smaller 

meteor counts at all 365 altitudes and the asymmetric meteor distributions seen in Figure 1 . RMS errors of these 

-1 

estimates 366 range from a few ms for the zonal components at central altitudes at SRO to -100% or greater 



367 uncertainties at YKF that render the latter useless in defining day-to-day variability in these 


368 quantities. 369 Estimates of semidiurnal tide zonal and meridional amplitudes for each day in Case 5 at 370 
3-km altitude intervals are displayed in the left and right panels of Figure 9. These estimates are 371 seen to be in 
close agreement with the specified values (dashed lines) at all but the lowest 372 altitude at TdF and the lowest and 

highest altitudes for KGI and BL. Estimates at SRO and YKF, 373 in contrast, are reasonably accurate at the central 

-1 

four altitudes, apart from sporadic departures of 374 ~20 to 60 ms at these altitudes, primarily at YKF and in the 
meridional component at SRO. 375 Hourly estimates of zonal and meridional GW momentum fluxes for the Case 5 
376 composite day at 3-km altitude intervals are shown in the left and right panels of Figure 10. 377 These are seen 
to be close approximations to the specified values at TdF for all but the lowest 378 and highest altitudes, and even to 
follow the semidiurnal momentum flux modulations at these 379 altitudes, but with larger uncertainties. The results 
for KGI and BL are comparable and also very 380 good, but exhibit somewhat larger fluctuations about the specified 
values than seen at TdF at the 381 central four altitudes altitudes and increasing uncertainties at KGI at 76.5, 79.5, 
and 97.5 km and 382 at BL at 76.5, 94.5, and 97.5 km. Composite day hourly momentum flux estimates for SRO and 
383 YKF are seen to occasionally follow the specified semidiurnal modulation at 85.5, 88.5, and 91.5 384 km, but 
with very large uncertainties. At higher and lower altitudes, these estimates exhibit very 385 large errors. 386 b.6. 
Case 6 387 We now compare the ability of the five radars to define mean and tidal motions and GW 388 momentum 
fluxes for the motion field defined by Case 6 in Table 2. This case includes mean and 389 tidal motions 
representative of higher latitudes and two transient GW packets having propagation 390 to the east and north and that 

occur randomly for 3 and 4 hr each day. The two GWs have zonal 391 and meridional propagation, daily mean 
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momentum fluxes of 50 and 25 ms , and periods of 20 392 and 30 min, respectively. Both also have horizontal and 
vertical phase variations with 393 wavelengths given in Table 2. Inferred monthly mean winds, tidal amplitudes, and 
GW 394 momentum fluxes obtained with each radar are shown in Figure 1 1 . Hourly estimates of zonal 395 and 
meridional GW momentum fluxes for a Case 6 composite day at 3-km altitude intervals are 396 shown in the left and 
right panels of Figure 12. Note here that the composite day momentum flux 397 variations are different for each 
radar in each component. Case 6 differs from Cases 1 and 4 in 398 having 1) different, but comparable, mean 



motions, 2) stochastic rather than uniform (and 399 stationary) GWs, and 3) vertical phase variations of the GWs. 

400 Monthly estimates of mean winds and tidal amplitudes shown in the three left columns of 401 Figure 1 1 have 
accuracies and biases that are very similar to those seen previously in Cases 1 402 and 4. Mean wind and diurnal tide 

amplitude estimates are again very accurate at TdF and nearly 403 as accurate at KGI and BL. Somewhat greater 

-1 

uncertainties are observed at SRO and YKF, but 404 even these are not larger than ~2 ms except at the lowest and 
highest altitudes. Semidiurnal tide 405 amplitude estimates are also nearly identical to those obtained in Cases 1 and 
4. Amplitudes are 406 again under-estimated by ~5% at KGI, SRO, and YKF, and by ~ 2 % at TdF and BL. The 
similar 407 accuracies achieved in these three cases suggest that the greater intermittency of the larger 408 amplitude 
GWs in Case 6 does not impact the ability of these radars to provide reasonably 409 accurate definitions of the 
monthly mean and tidal fields. 410 GW momentum flux estimates shown in the right column of Figure 1 1 are 
seen to be 41 1 relatively more accurate at KGI, TdF, and BL, where uncertainties are ~10% or less at central 412 
altitudes, but approach 20% at KGI and BL at the lowest and/or highest altitudes. Momentum 



413 flux estimates are also relatively accurate at SRO and YKF at the central 2 to 4 altitudes having 
414 the highest meteor counts, but exhibit large errors at lower and higher altitudes where meteor 415 counts are low. 
Accuracies at KGI, TdF, and BL are comparable to those seen in Case 1 for 416 which GW amplitudes and 
momentum fluxes are smaller. However, accuracies are not as good 417 as seen in Case 4 where the zonal 
momentum flux is largely determined by the stationary, large418 amplitude GW contributing the majority of the 
momentum flux. 4 1 9 Diurnal variability of the composite day zonal and meridional momentum fluxes for each 
420 radar shown in Figure 12 reinforce the statements above about the monthly mean profiles. The 421 TdF and BL 
radars capture the specified momentum flux distributions extremely well, but with 422 slightly greater uncertainties 
at TdF at the highest altitude, at TdF and BL at the lowest altitude, 423 and a greater loss of sensitivity at the highest 
altitudes at BL. Momentum flux estimates at KGI 424 are slightly less precise than either TdF or BL at central and 
lower altitudes, but are somewhat 425 better than at BL at the higher altitudes. In contrast, momentum flux estimates 
at SRO and YKF 426 are only accurate at the central ~4 and 2 altitudes, respectively, and are susceptible to very 
large 427 uncertainties below and above. These results suggest that high meteor counts (and more uniform 428 
sampling of the motion field) are essential in defining momentum fluxes when GWs are variable 429 and tides 
achieve large amplitudes. 430 In order to assess whether momentum flux estimates are negatively influenced by large 
43 1 tidal amplitudes, we also performed the same mean, tidal, and GW momentum flux assessments 432 for Case 6, 

but with zero mean winds and tidal amplitudes. This yielded 1) very accurate mean 433 wind estimates (except at the 
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lower altitudes at YKF), 2) diurnal and semidiurnal tide amplitude 434 estimates of ~1 ms at KGI, TdF, and BL and 

-1 

~2 ms at the central altitudes at SRO and YKF, 435 and 3) momentum flux estimates comparable to or slightly more 


accurate than with the tides 



436 present, with the greatest improvements at SRO. Hourly estimates for the composite day without 
437 mean winds and tides are likewise of comparable accuracy to those for the full field specified for 438 in Case 6 
in Table 2. 439 b. 7. Case 7 440 Case 7 is very similar to Case 6, with the same mean and tidal winds, but including 
two 441 intermittent large-amplitude GWs having propagation to the NE and SE. The GW propagating to 442 the NE 
occurs for 2 hr twice daily at random times; that propagating to the SE occurs for 1 hr 3 443 times daily at random 
times. The GWs have periods of 20 and 15 min and both horizontal and 444 vertical phase variations. Amplitudes 

and wavelengths of each are listed in Table 2. Together, 445 they result in mean zonal and meridional momentum 
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fluxes of 75 and -25 m s , respectively. 446 Inferred monthly mean winds, tidal amplitudes, and GW momentum 
fluxes obtained with each 447 radar for Case 7 are shown in Figure 13. Figure 14 shows hourly estimates of zonal 
and 448 meridional GW momentum fluxes for a Case 7 composite day at 3-km altitude intervals. 449 Monthly mean 
wind, tidal amplitude, and GW momentum flux estimates for Case 450 7 are very similar to those seen in Case 6 
above. Mean wind and diurnal tide amplitude estimates 451 are very accurate at KGI, TdF, and BL but exhibit larger 
errors at SRO and YKF at the highest 452 and lowest altitudes, respectively. Semidiurnal tides are under-estimated 
by ~5% at KGI, SRO, 453 and YKF, and by ~2% at TdF and BL, as in Case 6. However, larger errors are seen at 
SRO and 454 YKF at the highest and lowest altitudes, respectively, as seen for the mean winds and diurnal tide 455 
amplitudes. GW momentum flux estimates also exhibit comparable accuracies to those in Case 456 6, suggesting that 
GW packet duration does not have a strong influence on the ability to define 457 these quantities. Finally, as for Case 
6, we also performed these assessments with mean winds 458 and tidal amplitudes specified to be zero. Similar to the 

results for Case 6, this yielded 1) 459 extremely accurate mean wind estimates (except at the lower altitudes at YKF), 

-1 -1 

2) diurnal and 460 semidiurnal tide amplitude estimates of ~l-2 ms at KGI, TdF, and BL and ~2-5 ms at the 461 
central altitudes at SRO and YKF, and 3) momentum flux estimates comparable to or slightly 462 more accurate than 
with the tides present. As in Case 6, hourly estimates for the composite day 463 without mean winds and tides are of 
comparable accuracy to those for the full field specified for 464 in Case 7 in Table 2. 465 5 Summary and 
Discussion 466 We have assessed the relative measurement accuracies of meteor radars having similar 467 


frequencies (32.55 to 36.9 MHz) but different antenna configurations, peak power, and sampling 468 modes. Our 



evaluation was performed with meteor distributions observed by five radars 469 extending from 62.1 S to 62.5 N 
obtained during June of 2009, 2010, or 201 1 for which meteor 470 statistics tend to be better at northern than 
southern latitudes. Two of these radars, SAAMER on 471 Tierra del Fuego and DrAAMER on King George Island 

o 

(53.8 and 62.1 S, respectively) were 472 specifically designed for enhanced measurement capabilities for large-scale 
motions and GW 473 momentum fluxes. This was accomplished in each case with an 8-Yagi transmitting array 474 

o o 

directing power into 8 lobes centered at 35 zenith angles separated by 45 in azimuth. The other 475 three radars are 
conventional meteor radars at northern latitudes (see Figure 1). Our purpose was 476 to determine the relative 
abilities of different meteor radar configurations to recover specified 477 mean winds, diurnal and semidiurnal tide 
amplitudes, and GW momentum fluxes. 478 Seven test motion fields included various superpositions of constant 
mean winds, 479 constant and/or variable tidal amplitudes, and GWs having varying amplitudes, momentum 480 
fluxes, scales, periods, propagation directions, and intermittencies. The test fields were sampled 481 according to the 
observed meteor distributions throughout the month for each radar assuming no 482 radial velocity, range, or angular 

o 

uncertainties in the measurements. Only meteors at zenith 483 angles between 15 and 60 were used, due to the large 
altitude uncertainties accompanying ~1 484 zenith angle uncertainties in angle-of-arrival estimates at larger zenith 
angles. The various fields 485 were estimated in 3-km altitude bins (centered from 76.5 to 94.5 km) using the method 
described 486 by FI 0b, in which mean winds and tidal amplitudes are determined by S-transform fits and 487 
removed from each radial velocity distribution before application of the Hocking [2005] method. 488 For each case 
mean wind, tidal amplitude, and GW momentum flux profiles were determined for 489 monthly and/or daily 
intervals (as appropriate). Also estimated in cases having diurnal variations 490 in GW momentum fluxes were 
hourly profiles throughout a composite day. 491 Our evaluation of relative radar performance revealed the following: 
492 1) measurement accuracies depend strongly on meteor counts and antenna beam patterns; 493 accuracies 
improve with higher meteor counts, smaller zenith angles, and more 494 symmetric beam patterns enabled by 

o 

crossed antennas, 495 2) SAAMER (on TdF), having -12,000 meteors/day between 15 and 60 zenith angles 496 


(-120 meteors/hr in a 3-km altitude bin at -90 km) performed best overall, yielding quite 497 accurate mean wind, 



tidal amplitude, and GW momentum flux estimates for all test cases, 498 3) the BL meteor radar and DrAAMER (on 
KGI) performed comparably and very well, 499 despite DrAAMER having -30% smaller meteor counts (-8,000 
meteors/day compared 500 to -12,000 total at BL, or -80 and 120 meteors/hr in a 3-km altitude bin near the peak of 
501 the meteor distribution); however, DrAAMER yielded a consistently larger under502 estimate of semidiurnal 
tide amplitudes (-5% rather than the -2% seen at BL), and 503 4) the SRO and YKF meteor radars having 6 kW 
peak power and antenna with only one 504 polarization (having - 5,000 and 3,000 meteors/day, respectively, or 50 
and 30 505 meteors/hr in 3 km) provided less accurate measurements and over more limited 506 altitudes; they 
nevertheless described monthly mean winds and tidal amplitudes 507 adequately, with the best results at central 
altitudes; monthly mean and hourly composite508 day GW momentum flux estimates were generally poor at lower 
and higher altitudes, but 509 reasonable at -85 to 90 km. 510 More general conclusions obtained from 
inter-comparisons of our various cases include 5 1 1 the following: 512 5) momentum flux estimates are more 

accurate when tidal amplitudes are small (or zero) 513 and when GW amplitudes and momentum fluxes are large, 
514 6) momentum flux estimates are not significantly impacted by complexity of the GW field, 515 including 
superposition, multiple scales and frequencies, and intermittency, 516 7) well configured radars (those having 
symmetric beam patterns ( e.g crossed Yagis) 517 and/or more meteors at smaller zenith angles) capture daily and 
composite-day hourly 518 variability of total winds and GW MFs well at altitudes where meteor count are high, and 
519 8) radars that have asymmetric beam patterns (crossed Yagis) and/or low meteor counts 520 (less than -25 
meteors/hr in 3 km) will likely exhibit significant biases in estimates of 521 daily mean winds and tidal amplitudes; 
they will likely also exhibit significant biases in 522 estimates of hourly composite-day GW momentum fluxes and 
variability of these fields 523 on longer time scales. 524 Our assessment of meteor radar performance above 
assumed no radial velocity, range, or 525 angle-of-arrival measurement uncertainties. Introducing such uncertainties 


would of course 526 degrade measurement accuracies to some degree in all cases. Here we assess the likely impacts 



527 of such measurement uncertainties for the magnitudes of the mean, tidal, and GW motions and 
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528 momentum fluxes considered above. 529 If we assume that RMS radial velocity uncertainties are ~1 to 2 ms 
(which is a 530 reasonable estimate for SAAMER), then these contributions to measured radial velocities are 531 
significantly less than typical radial velocity magnitudes, given the mean winds and tidal and 532 GW amplitudes 

o 

employed for our assessment. For meteors at the largest zenith angles, -60 ,533 these uncertainties would contribute 
little to GW momentum flux assessments given their small 534 implied errors in vertical velocity contributions to 
radial velocity variance relative to the large 535 contributions by large horizontal velocities. These uncertainties at 
large zenith angles would also 536 contribute very small errors in estimates of mean and tidal motions, which are 
obtained by 537 averaging large horizontal motions rather than effectively differencing radial velocity variances 538 
for GW momentum flux assessments. 539 At smaller zenith angles, there is a greater potential for influences of radial 
velocity 540 uncertainties, given that these angles provide greater sensitivity to the differential radial 541 velocities 
with azimuth that contribute to the momentum flux determination. The most sensitive 542 radars to such 
uncertainties will be those focusing transmitted power at small zenith angles, 543 specifically SAAMER and 

o 

DrAAMER, both of which achieve maximum meteor counts at ~35 544 zenith angles. 545 To assess qualitatively 
the impacts of radial velocity uncertainties on measurements of 546 GW momentum fluxes, we assume that all 

o 

meteors occur at the 35 zenith angle of maximum 547 meteor counts for SAAMER and DrAAMER. We then 
evaluate the relative uncertainty of the 548 momentum flux assessment using the “dual-beam” analysis of Vincent 


and Reid [1983] for a GW 



549 propagating zonally with superposed tidal winds and nominal radial velocity measurement 
550 uncertainties. The east and west radial velocities are given by 55 1 Ve = (u’ + Ud + Us)sin0 + w’cosO Ave’ 

(4) 552 and 553 Vw = -(u’ + Ud + Us)sin0 + w’cosO Avw’ (5) 554 where u’ and w’ are the horizontal and vertical 
GW velocities, Ud and Us are the zonal diurnal 555 and semidiurnal tidal motions (assuming negligible vertical tidal 

o 

motions and no horizontal 556 phase variations), 0 = 35, Ave’ and Avw’ are the radial velocity uncertainties, and 
angle 557 brackets denote temporal averaging. Assuming all of the velocities in Eqs. (4) and (5) are 558 uncorrelated 

except for u’ and w’, which are in phase or antiphase, squaring, subtracting Eq. (5) 559 from Eq. (4), and 
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rearranging, the measured momentum flux may be written as 560 <u’w’>m = (Ve - Vw + Ave’ -Avw’ )/2sin(20) 

2 2 

(6) 561 = <u’w’>t + (Ave’ -Avw’ )/2sin(20) 562 where <u’w’>t is the true GW zonal momentum flux for the 

specified motion field. The 563 fractional uncertainty in the measured momentum flux is then 564 A<u’w’>/<u’w’>t 

22 

= (Ave’ -Avw’ )/[<u’w’>t 2sin(20)] (7) 565 Thus, fractional uncertainties are minimized if radial velocity 
uncertainty variances are small 566 relative to the true GW momentum flux, as 2sin(20) = 1.88 for the zenith angle 

specified. The 567 largest magnitude of the numerator occurs if the radial velocity uncertainty variance is negligible 
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568 in one beam. But for the uncertainties assumed above, this is Av’ -3 m s , with likely values 569 significantly 
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smaller. By comparison, mean GW momentum fluxes in the MLT are typically -10570 ms , with peak values -3 to 
10 times larger, at most locations, with values over the Drake 571 Passage “hotspot” expected (and measured, see 


FI 0b) to be -2 to 5 times larger. These values 



572 suggest a maximum fractional uncertainty of the momentum flux due to radial velocity 
573 uncertainties of -15% at typical sites and ~5% in GW source hotspots. 574 We anticipate range uncertainties 

o o 

comparable to the range resolution of 2 km and angle575 of-arrival uncertainties of - 1 in zenith angle and ~1.5 in 
azimuth. These imply that altitude 576 uncertainties are defined by zenith angle uncertainties near zenith angles of 

o 

-60 and by range 577 uncertainties at smaller zenith angles. At all zenith angles, however, these values suggest 
altitude 578 uncertainties comparable to or smaller than our chosen 3 -km altitude bin. Thus range and angle579 
of-arrival uncertainties appear unlikely to significantly impact measurements of mean and tidal 580 winds having 
large vertical scales or GW momentum fluxes that rely most on radial velocities 581 occurring at smaller zenith 
angles. These results suggest that measurement uncertainties are 582 likely to contribute much less to momentum 
flux uncertainties than inadequate definition of the 583 mean, tidal, and GW fields due to low meteor counts arising 
due to low radar sensitivity or small 584 averaging intervals and/or range bins. 585 Several previous studies have 
assessed the potential for momentum flux measurement 586 uncertainties due to inadequate sampling of the motion 
field. Kudeki and Franke [1998] and 587 Thorsen et al. [2000] evaluated momentum flux measurements employing 
the “dual-beam” 588 technique and concluded that very long averaging intervals were required to achieve statistical 
589 significance. Kudeki and Franke [1998] assumed the GW velocity fields to be defined by 590 Gaussian 
distributions with a net momentum flux equal to 1% of that assuming a single GW was 591 present and inferred a 
required integration time of 16 days. This assumption, however, is in 592 contradiction to the strong correlations 
among component velocities that often occur when one or 593 several large-amplitude GWs dominate the 
high-frequency motion field and the total momentum 594 flux, as is often observed. In such cases, the averaging time 


implied by the Kudeki and Franke 



595 [1998] would decrease by 2 decades or more. Examples of dual-beam or multiple-beam 

596 measurements employing MF, VHF, and UHF radars that provide clear evidence for significant 597 and variable 
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GW momentum fluxes (magnitudes as large as ~70 ms) occurring on time scales 598 as short as a few hours or less 

include the following: 599 1) MLT measurements with the MF radar at Buckland Park, Australia showing clear 
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diurnal 600 tide modulation of GW momentum fluxes with a modulation amplitude of ~30 m s 601 [Fritts and 

Vincent , 1987], 602 2) MLT measurements with the SOUSY VHF radar at Andoya, Norway revealing 603 
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momentum fluxes (per unit density) of ~66 m s over ~3 hours [Reid et al., 1988], 604 3) troposphere and 
stratosphere measurements using the MU VHF radar at Shigaraki, Japan 605 for multi-beam measurements of GW 
momentum fluxes over 6 days that exhibit clear 606 episodic enhancements as short as an hour simultaneously in 
multiple beams and 607 spanning multiple altitudes [Fritts et al., 1990], 608 4) MLT measurements employing the 
former VHF radar at Poker Flat, Alaska that revealed 609 clear anti-correlations between diurnal and semidiurnal 

tide winds and GW momentum 610 fluxes over 4 and 1 6 day intervals exhibiting stronger momentum flux 
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modulations (~5 to 61 1 10ms) accompanying larger tidal amplitudes [ Wang and Fritts, 1991], 612 5) MLT 

measurements spanning two 10-day intervals with the Jicamarca VHF radar in Peru 613 that exhibit daily mean 
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momentum fluxes of ~10 m s and maxima of ~30 m s , with 614 hourly profiles exhibiting significant coherence in 
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altitude and time and maximum 615 magnitudes of ~60 m s [ Fritts et al., 1992], 



616 6) MLT measurements with the MF radar at Buckland Park, Australia exhibiting clear 


617 diurnal modulation of GW momentum fluxes by a large-amplitude 2-day wave [ Murphy 618 and Vincent, 1998], 

619 7) MLT measurements employing the dual-beam Arecibo UHF radar that revealed variable 620 GW momentum 
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fluxes often near zero, but exhibiting occasional maxima as large as -50 621 m s , that were largely anti-correlated 
with the large-scale wind field [Fritts et al., 622 2006b], and 623 8) MLT measurements spanning ~8 hr employing 
the Arecibo UHF radar that revealed very 624 significant GW activity, but with momentum fluxes that were very 
small, suggesting 625 ducting [Fritts and Janches, 2008]. 626 A more recent study by Vincent et al. [2010] 
employed numerical simulations of test GW 627 fields with a statistical model of radial velocity and angle-of-arrival 
uncertainties to evaluate the 628 accuracy of meteor radar measurements of mean winds and GW velocity variances 
and 629 momentum fluxes as functions of the meteor rate within a 2-km range bin. These authors 630 concluded that 
mean winds could be determined with relatively few meteors, but that estimation 63 1 of GW variances, and 
especially momentum fluxes, with small uncertainties required 632 considerable averaging enabling large meteor 
counts. The meteor counts employed to assess 633 measurement capabilities in the study ranged from 10 to 200. 634 
By comparison, our assessment above assumed a need for significant meteor counts, with 635 the typical averaging 
interval for definition of mean and tidal winds and GW momentum fluxes 636 being a month, but with mean wind 
and semidiurnal tide assessments also evaluated for one day 637 or for one hour of a composite day for the month 
(hence 24 or 30 hr of data). Employing real 638 meteor distributions observed during a test month for each of five 


radars ( June of 2009, 2010, or 



639 2011), these yielded meteor counts in a 3-km altitude bin at the peak of the meteor distribution 
640 for the month of -80,000 at TdF, -70,000 at BL, -40,000 at KGI, and -25,000 at SRO and YKF. 641 Meteor 
counts at the upper and lower edges of these distributions varied by radar, but were 642 typically smaller by -2 to 4 
times. Daily and composite-day hourly assessments employed -3 643 and 4% of these meteor counts, respectively. 
Thus, in all cases, our meteor counts for each radar 644 were significantly larger than those assessed by Vincent et al. 

[2010] near the peak of the meteor 645 distribution. 646 Assuming that accuracies vary as the signal-to-noise ratio 

1/2 

(S/N) and thus as N , where N 647 is the meteor count, we can compare our inferred accuracies with those of 
Vincent et al. [2010], 648 As noted above, the radars at BL, TdF, and KGI have peak monthly meteor counts of N 
-40,000 649 to 80,000 in a 3-km altitude bin, with N -10,000 to 25,000 at the edges of the distributions. 650 These 

imply reductions by factors of -15 to 25 in the mean wind measurement errors at the 651 meteor distribution peak 

-1 

relative to the 5 ms uncertainties displayed by Vincent et al. [2010] in 652 Table 1, with reductions by factors of -10 
at the highest and lowest altitudes. Inspection of the 653 monthly mean winds for the seven cases considered above 

suggests that our results are largely 654 consistent with that expectation. Similarly, we expect error reductions of -4 

-1 

and 2 in the center 655 and edge regions for daily mean winds, implying errors of -1 to 3 ms . These are largely 656 

consistent with the daily mean winds inferred for our Case 5 test fields displayed in Figure 8. 657 Finally, meteor 

-1 

counts at SRO and YKF suggest monthly mean wind uncertainties smaller than 658 the -5 ms uncertainties of 
Vincent et al. [2010] by -12 and 5 in the center and edge regions. In 659 these cases, our uncertainties seem to be 
comparable to, or somewhat larger than, those expected. 660 We note, however, that the antenna patterns implied by 
the meteor locations shown in Figure 1 661 suggest significant asymmetries that likely also contribute to 


measurement errors. 



Applying the same factors to our momentum flux estimates, we expect to see errors of 


663 ~10 and 25% for the monthly means at the center and the edges of these distributions. Again, our 664 results are 
consistent with these expectations, though with variations that depend on the 665 complexity of our motions fields 
and the meteor distribution of each radar. Note, in particular, 666 that Vincent et al. [2010] did not include tidal 
winds in their assessment, but that variable tidal 667 winds can induce apparent mean winds and GW momentum 
fluxes if the meteor distribution is 668 not approximately uniform in space and time. Our results also are relatively 
more accurate for 669 cases having constant GW character throughout the month; they are less accurate for Cases 6 
and 670 7 for which GWs are randomly distributed and likely more similar to the GW spectra specified 671 by 
Vincent et al. [2010], We cannot assess the contributions of the assumed radial velocity and 672 angle-of-arrival 
uncertainties employed by Vincent et al. [2010] to the overall uncertainties 673 displayed in their Table 1. However, 
our own assessment of radial velocity uncertainties above 674 suggests that the majority of the measurement errors 
seen in our results can be attributed 675 primarily to the expected dependence on meteor counts. 676 677 6 
Conclusions 678 Our assessment of meteor radar measurement capabilities presented in previous sections 679 has 
demonstrated both 680 1) clear capabilities for quantitative measurements of mean winds, tides, and GW 681 
momentum fluxes by radars that achieve sufficiently high meteor counts and provide 682 sensitivity to radial 
velocities at sufficiently small zenith angles, and 683 2) strong dependence of these measurement capabilities on 
radar power, beam geometry and 684 sensitivity, and spatial and temporal meteor distributions. 



685 Meteor radars that combine high power, high meteor counts at small zenith angles, and 


686 symmetric beam patterns exhibit the smallest measurement errors. Those having low power, 687 acquiring a 
majority of meteor detections at large zenith angles, and/or experiencing beam or 688 meteor detection asymmetries 
exhibit much larger measurement uncertainties, thus either 689 requiring longer averaging intervals or precluding the 

more challenging measurements, 690 especially estimations of GW momentum fluxes. 691 Monthly mean winds 

-1 

can typically be measured to within 1 ms over ~15 to 20 km, with 692 maximum accuracies between ~85 and 90 

km, for radar frequencies near 35 MHz and having 693 higher power that conventional meteor radars. These radars 

-1 

can also provide daily mean winds 694 having accuracies of ~1 to 3 ms where meteor counts are sufficient, even in 
the presence of 695 large tidal motions. Monthly mean diurnal and semidiurnal tide amplitudes are also recovered 
696 very well by higher-power meteor radars, but typically with amplitude under-estimates of a few 697 %. Mean 
winds and tidal amplitudes are recovered much more accurately, however, by meteor 698 radars that achieve high 
meteor counts at small zenith angles compared to meteor radars 699 achieving the same meteor counts, but at larger 
zenith angles. 700 GW momentum flux estimates are challenging for meteor radars because of the large 701 
number of meteors required to adequately resolve differences in radial velocity variances at 702 opposite azimuths. 
Nevertheless, accurate estimates over sufficiently long time scales (typically 703 averages over ~1 to 30 days) appear 
to be possible with sufficient meteor counts (~2,000 to 704 60,000 in an altitude bin). Accuracies are also enhanced 
with high meteor counts at small zenith 705 angles, given the greater sensitivity of these radial velocities to the 
correlations between 706 horizontal and vertical velocities within the GW field. Importantly, meteor radars appear to 
707 measure momentum fluxes for stationary and propagating GWs, for complex GW 708 superpositions, and for 
GW fields exhibiting significant intermittency equally well on monthly 709 time scales. 710 Acknowledgments 
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Table 1. Radar characteristics for the five radars employed in this evaluation of measurement 880 capabilities. 
Meteor counts for zenith angles between 15 and 50 are ~50% of the total counts 881 (right column) for KGI and TdF 


and ~30% of the total counts for the other radars. 


882 883 Table 2. Mean, tidal, and GW parameters used for test cases evaluating radar measurement 884 capabilities 




employing real meteor distributions and test motion fields. Mean GW momentum 
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Figure 1. 
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illustratin 
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for the 
five radars 
employed 
908 in this 
assessmen 
t of 

measurem 

ent 

capabilitie 
s. The 
radar 
locations 
and 

frequencie 
s are listed 
in 909 
each 

panel. The 
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BL, TdF, 
and KGI 


use 



crossed 
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and have 
nearly 910 
symmetric 
beam 
patterns; 
The radars 
at YKF 
and SRO, 
however, 
use 
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(linearly- 

911 

polarized) 
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h and 

east-west, 

respective 
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resulting 



in the 
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and 
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measurem 

ents, 

respective 
ly, at 
those 913 
sites. 

Nulls at 
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are due to 
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of meteor 



detections 


that 914 
may have 
ground 
clutter 
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tion. 915 
Figure 2. 
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employed 
for mean 
and tidal 
wind and 
momentu 
m flux 
estimates 
916 for 
each radar 
(labeled at 
left). 
Meteor 
distributio 
ns are 


shown as 



(left) 
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counts in 
each 917 
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range bin, 
(center) 
counts per 
day, and 
(right) 
counts per 
hour for a 
composite 
day. 

These 918 
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angles 
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15 and 60 
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examined 



here. 919 


Figure 3. 

Monthly 

mean 

(left) 

winds, 

(second 

column) 

diurnal 

tide 

amplitude 
s, (third 
column) 
920 

semidiurn 
al tide 
amplitude 
s, and 
(right) 
GW 

momentu 
m fluxes 
for Case 
1. Solid 


and 



dashed 


921 lines 
denote 
zonal and 
meridiona 
1 profiles 
and the 
radars are 
designated 
at the left 
of each 
row. 922 
Vertical 
dashed 
lines show 
specified 
values. 
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Figure 4. 

As in 
Figure 3, 
but for 
Case 2. 
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Figure 5. 



As in 


Figure 3, 
but for 
Case 3. 
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Figure 6. 

As in 
Figure 3, 
but for 
Case 4. 

926 

Figure 7. 

As in 
Figure 3, 
but for 


Case 5. 



927 Figure 8. Daily (left) zonal mean winds, (second column) meridional mean winds, (third 928 column) zonal 
diurnal tide amplitudes, and (right) meridional diurnal tide amplitudes for Case 5. 929 Vertical dashed lines show 
specified values. 930 Figure 9. Daily (left) zonal and (right) meridional semidiurnal tide amplitudes for Case 5 for 
the 93 1 five radars (labeled at the top of the left panels). Colored lines show measured amplitudes from 932 (black) 
76.5 to (red) 97.5 km, dashed lines show specified amplitudes. 933 Figure 10. Hourly (left) zonal and (right) 
meridional GW momentum fluxes for Case 5 for a 934 composite day for the five radars (labeled as in Figure 9). 
Colored lines show measured 935 amplitudes from (black) 76.5 to (red) 97.5 km. Dashed lines show specified 
amplitudes. 936 Figure 11. As in Figure 3 for Case 6. 937 Figure 12. As in Figure 10, but for Case 6. 938 Figure 
13. As in Figure 3 for Case 7. 939 Figure 14. As in Figure 10, but for Case 7. 
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